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1. Introduction 

At high temperatures strongly interacting matter undergoes an transition to a new state where 
quark and gluons are not subject to confinement; it shows up in a qualitative change in the equation 
of state (see [jj]] for a review). The most common way to calculate the QCD equation of state is 
the integral method. The pressure is calculated as an integral of the trace of the energy momentum 
tensor e — 3p or the trace anomaly 
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The study of the trace anomaly is also interesting in its own right. At high temperature it is de- 
termined by the running of the QCD coupling constant a s , and at leading order it is proportional 
a s 2 T 4 . Thus it provides a sensitive test of the weakly coupled nature of the quark gluon plasma. For 
this reason it is also often called the interaction measure. Thermodynamic quantities are expected 
to have large cutoff effects at high temperatures when calculated using an unimproved fermion 
formulation, even if the lattice spacing is small in absolute units. These cutoff effects arise from 
the distortion of the quark dispersion relation on the lattice and can be cured by using an improved 
action [Q, g]. On the other hand, at low and intermediate temperatures the cutoff effects in ther- 
modynamic quantities can be understood in terms of the cutoff effects of the hadron spectrum. In 
particular, the breaking of the so-called taste symmetry in the staggered fermion formulation could 
be the dominant source of cutoff effects when using staggered quarks. To control cutoff effects in 
pressure and other thermodynamic quantities, one first needs to understand the cutoff dependence 
of the trace anomaly. In this contribution we will study in detail the cutoff dependence of the trace 
anomaly and compare our lattice results with the hadron resonance gas model at low temperature 
and with resummed perturbation theory at high temperatures. 
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Figure 1: The difference in temperature (left) set by the r\ and fK scales, described in the text. The running 
of the gauge coupling (right). 



2. Numerical results 



We calculate the interaction measure in 2+1 flavor QCD using a tree-level improved action 
for gauge fields and the highly improved staggered quark (HISQ) action [Eh. This combination 
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of gauge and quark action is called the HISQ/tree action. It can largely reduce cutoff effects in 
thermodynamic quantities both at low and high temperatures I^q]. 

The trace anomaly can be written in terms of the expectation values of gauge action and quark 
condensates as 

r4 = Rp [{S g auge)o ~ {S gauge ) T } - RpR m [2mi((ll} - (ll) T ) + m s ({ss)o - (ss) T )], (2.1) 

where the subscript "0" refers to T = and "T" to finite temperature for observables evaluated at the 
same value of the cutoff. Furthermore, Rp and R m are the lattice beta function and mass anomalous 
dimension that are defined below. The above expression is free of ultraviolet divergences. On a 
hypercubic lattice x N T the physical temperature is set by the size of the temporal dimension 
and the lattice spacing as T = 1/ (N T a). For T = calculations we use N T ^ A^, and for T > we 
keep N S /N T = 4 and at fixed N T vary the lattice spacing a by varying the gauge coupling j8 = 10/ g 2 . 
The continuum limit in this setup corresponds to N t —> °°. Therefore, we carried out this study on 
lattices with N t = 4, 6, 8, 10 and 12. The strange quark mass m„ is tuned to the physical value, 
while the two degenerate light quarks have masses mi = m s /20, slightly heavier than physical 
(mi ~ m v /27). In the continuum limit these light quark masses correspond to the pion mass of 
about 160 MeV. To set the lattice spacing in physical units (fm) we use the r\ scale defined through 
static quark potential 

r Z ^\ r=n = h (2.2) 

We use the value r\ = 0.3106 fm [Q]. Alternatively, the kaon decay constant, /«• = 156.1 MeV is 
used to set the scale. >From Fig. [j] (left) one can see how the choice of reference observables affects 
the conversion of temperature from lattice units to MeV. Over the temperature range of interest on 
N r = 6 lattices the difference is within 9 MeV, and on N t = 12 within 2 MeV. Another effect of 
using different scales is related to the /3 -function and mass anomalous dimension that enter in Eq. 

where m s (j8) defines a line of constant physics (LCP), i.e., the combination of the gauge coupling 
and strange quark mass such that the kaon mass (in MeV) stays approximately constant in the 
whole j8 range used in the simulation. Rp for r\ and fx scale is shown in Fig. [j] (right), together 
with the 2-loop perturbative result. 

The present status of the interaction measure with the HISQ/tree action is shown in Fig. || 
together with the stout continuum estimate of Ref [g]. For temperature T < 180MeV we see good 
agreement between our HISQ/tree results and the stout continuum estimate. At higher tempera- 
tures we see significant discrepancies between HISQ/tree results and the stout results, though the 
discrepancies seem to disappear at temperatures T > 350MeV. The height of the peak decreases 
slightly when going from N t = 8 to N t = 10 and 12. Thus one may expect that the discrepancies 
between HISQ/tree and stout actions in the peak region will be reduced once the continuum extrap- 
olation for HISQ/tree is performed. It should be noted that in the previous calculations with p4 and 
asqtad actions [|9[ [R], [TTJ] the trace anomaly was significantly smaller due to large cutoff effects. 

Next we study the effect of the scale setting on the trace anomaly. In Fig. |2| (right) the results 
for N t = 6 and 10 are shown using r\ and to set the scale. Different scale settings affect the 
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Figure 2: Comparison of the HISQ/tree interaction measure on N T = 4, 6, 8, 10 and 12 lattices with the 
stout continuum estimate [||] (left), N? = 6 and 10 HISQ/tree data (right). Filled (open) symbols in the right 
panel correspond to the r\ (fx) scale, see text. 

N T = 6 results, in particular in the peak region the trace anomaly is smaller if fx is used to set the 
scale. However, for N z = 10 there is almost no difference between the two scale setting procedures. 

We also calculated the trace anomaly using the asqtad action. The expression for the trace 
anomaly in terms of expectation values of local operators is given in Ref. JTOp- In Fig. [3| we 
compare our results obtained with HISQ/tree action for N t = 8 and N t = 12 using the scale set 
with both r\ (filled symbols) and fx (open symbols). If r\ is used to set the scale, there is a 
significant discrepancy between the HISQ/tree and the asqtad result for T < 200 MeV for /V T = 8. 
As discussed above this is due to large cutoff effects at low temperatures when the asqtad action 
is used. The discrepancy is significantly smaller for = 12, as expected. If we use fx to set the 
scale, the discrepancies between HISQ/tree results and asqtad results are greatly reduced. Note, 
however, that the peak height in the N T = 8 asqtad data is reduced when fx is used to set the scale. 
Fig. H also shows that the difference between the two scale setting procedures for the asqtad action 
is reduced in the case of /V T = 12 data. This is again expected. There is, however still a significant 
reduction in the peak height if fx scale is used. 

The low-temperature behavior of the interaction measure is shown in Fig. || (left) and com- 
pared to the hadron resonance gas (HRG) calculation shown as the black line. At low temperatures, 
cutoff effects arising from the taste symmetry breaking of the staggered fermion formulation could 
be significant. To take into account possible cutoff effects in the comparison with HRG we replaced 
the contribution of each pion and kaon with a contribution averaged over the sixteen different tastes 
of pseudo-scalar mesons. The pseudo-scalar meson masses for each taste as function of the lattice 
spacing have been previously estimated by the HotQCD collaboration [g]. In Fig. || we show the 
HRG results with the modified pseudo-scalar meson sector as colored lines for each value of N T . 
We use the same color coding for the lattice results and the modified HRG result. Despite the fact 
that even for HISQ/tree action the effects of taste breaking in the pseudo-scalar sector are signifi- 
cant, the resulting cutoff effects are surprisingly small. They are of the same order or smaller than 
the statistical errors for all N T values, including N T = 6. Fig. ^| shows that the lattice result starts to 
disagree with the HRG model at T ~ 150 - 160 MeV. 

Finally, let us compare our results for the trace anomaly with earlier calculations using the p4 
action H [RJ as well as with resummed perturbation theory [12]. In Fig. ^| (right) we compare 
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Figure 3: The interaction measure for the HISQ/tree action on N T = 8 (top) and N T = 12 (bottom) ensembles 
with r\ (left) and fx (right) scale. 

our results for the HISQ/tree action with previous results obtained with p4 action and the stout 
continuum estimate. For T > 400 MeV different lattice data agree well with each other. The only 
exception is the N t = 4 p4 data, which show the expected cutoff effects at high temperatures. The 
N t = 4 HISQ/tree data should show similar cutoff dependence, but interestingly enough this is not 
observed. It is important to stress that at high temperatures the cutoff effects for p4 action should 
be small for N t > 6. Thus the agreement between p4 results and HISQ/tree results is expected. We 



see quite good agreement with the resummed perturbative results [12] 



3. Conclusion 

We studied the trace anomaly in QCD using HISQ/tree and asqtad actions and lattices with 
temporal extent A^ T =4, ,6, 8, 10 and 12. Using different observables to set the scale, r\ and fx, 
allows for a crude estimate of the magnitude of cutoff effects, which seem to be small at the finest 
/V T = 12 lattices both for HISQ/tree and asqtad actions. Clearly to get reliable continuum results for 
the trace anomaly and thus for the equation of state the calculations on Af T = 10 and 12 lattices need 
to be significantly extended. We also compared our lattice results with HRG at low temperatures 
and resummed perturbative results at high temperatures. The trace anomaly is well described by 
HRG for T < 150MeV. We also find a good agreement between lattice and resummed perturbative 
results at high temperatures. 
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Figure 4: The interaction measure at low (left) and high (right) temperature. The curves on the left panel 
represent the HRG calculation, see text. On the right the dashed curve shows the resummed perturbative 
result with 2-loop running and the solid curve the resummed perturbative result with 1-loop running. 
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